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Large twist-angle grain boundaries in layered structures are often described by 
Scherk's first surface whereas small twist-angle grain boundaries are usually described 
in terms of an array of screw dislocations. We show that there is no essential distinc- 
tion between these two descriptions and that, in particular, their comparative energetics 
depends crucially on the core structure of their screw-dislocation topological defects. 
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Contrary to their name, minimal surfaces are widely observed. It has been proposed 
that lamellae in diblock copolymer systems form minimal surfaces so as to minimize the 
interfacial area between blocks. Other systems, such as cubic phases @ and the smectic-Q 
phase are also described in this way. These minimal surfaces contain grain boundaries, 
which rotate the layer normals. On the other hand, a twist-grain-boundary (TGB) in a 
smectic-A liquid crystal phase is described in terms of a lattice of screw dislocations |3| . We 
will show that these two different descriptions lead to nearly identical geometrical surfaces 
even though the philosophy behind their construction is quite different. 

Minimal surfaces are simply surfaces of zero mean curvature. Because of this, they 
play an important role in the physics of complex fluids: their curvature energy always 
vanishes. Though they may have periodically repeated unit cells or, like Scherk's first 
surface f|, provide a continuous transition between parallel planes at different angles, 
they do not have an inherent length scale - a uniform dilation of a minimal surface is 
still minimal. One of the simplest minimal surfaces is the helicoid: a surface with height 
function h(x,y) = h + [6/(27r)j tan _1 (y/x). This surface is smooth and continuous, and 
there are no constraints on the parameter b giving the pitch of the helicoid. A screw 
dislocation in a smectic-A, however, is a topological defect in a lamellar structure with 
a specified periodic spacing d. It has a height function identical to that of a helicoid, 
though its pitch must be an integer multiple of the layer spacing. As we discuss below, 
this constraint is most easily implemented in an Eulerian picture in which the single- valued 
phase of the mass-density wave is spatially dependent. Moreover, a true dislocation has a 
singular core because the layer spacing and thus the strain diverges there. This costs energy 
precisely because there is a preferred layer spacing. Thus the existence of a well-defined 
periodic reference state distinguishes screw dislocations from helicoids. In this letter, we 
will explore the consequences of this distinction. 

We start by considering Scherk's first surface. The height function h(x, y) is given 
implicitly via ||: 

cos [cos (%) h(x, y) + sin (a/2) y\ 

x sin a = In p -. — r — — — — \ 

cos [cos (^j h(x, y) — sin (ct/2) y\ 

where t = (x, y, h(x, y)) and Ni^ = (0, ± sin(a/2), cos(a/2) are unit vectors. By using 
the rightmost expression in (|l|), it is easy to see that as x — ► — oo (for < a < ^) the 
surface is the solution of cos (Ni ■ t) = 0, i.e., Scherk's surface is a family of parallel planes 
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with layer normal Ni and spacing d = 7r. Similarly, as x — > oo Scherk's surface becomes a 
family of parallel layers with normal N2 and spacing d = n. 
Note that as a — * 0, ([[]) becomes 



which is a single helicoid with pitch p = 2it = 2d located at the origin. Since planes 
at infinity are separated by tt, this can be interpreted as a +2 screw dislocation. Note, 
however, that the normal distance between the planes diverges along the fact that 

will prove to be essential in our discussion of energetics. In biphasic materials, such as 
diblock copolymers, these layers should be alternately identified as A-blocks and B-blocks. 
This would make the layer spacing at infinity d = 2d = 2n and thus the Burgers vector 
b = 2n would be a strength +1 dislocation However, in the following we will take each 
layer to be identical. It will only be convention as to the strength of the dislocation - the 
actual length of the Burgers vector is the appropriate invariant. 

For finite a, the intersection of the Scherk's surface with the plane x = consists 
of two sets of equally spaced lines: one parallel to the 5-axis at positions yk = kid for 
integer k, with separation Id = tt/ sin(a/2) = dj sin(ct/2) and a second set parallel to 
y-axis at the heights hk — kd' with separation d' = %/ cos(a/2) = dj cos(a/2). Around 
the first set of lines, Scherk's surface is a local helicoid with height function h(x, y) = 
— [1/ cos(a/2)] tan -1 [x cos(a/2)/5y] for small 5y = y — yk- Thus these helicoids have 
pitch p z = 2d' equal to twice the separation between the planes parallel to y at x = 
0. Similarly, around the second set of lines there are helicoids: for small Sh = h — hk, 
y(x,h) = — [1/ sin(a/2)] tan -1 [xsm(a/2)/5h]. Thus they have pitch p y = 2ld equal to 
twice the separation between planes parallel to z at x = 0. Note that neither p z or p y 
are integer multiples of the layer spacing at infinity d. The two sets of helicoids can be 
interchanged through the "dual" transformation h(x,y) y and a —>■ n — a which leaves 
Scherk's surface invariant and interchanges d' and Id- In the following we shall show that 
if the (non-chiral) ground state of the system is lamellar, with a preferred layer spacing, 
then energetics will break the geometric degeneracy of the dual mapping by converting one 



x sin a « In [1 — 2y sin (a/2) tan/i] w —y sin a tan h 



(2) 



so that 




(3) 



2 



set of helicoids into true topological screw dislocations with cores. Finally, for comparison 
with the smectic structures we will consider, we calculate the derivatives of h(x,y): 

d x h(x,y) =~ sin (a/2) ^^-LpO) (4) 
d y h(x, y) = tan (a/2) ^gff^afl) 

where 9 = y sin (a/2) and 7 = \xsirva. From these derivatives and boundary conditions 
Scherk's surface can be reconstructed. 

We will now consider a surface constructed via a linear superposition of screw disloca- 
tions (LSD) in a smectic-A phase, parallel to z with separation la along the y-axis. Recall 
that a smectic is described in terms of a complex scalar mass-density wave 

where q = 2ir/d and d is the layer spacing. While the surface h(x,y) is a true height 
function expressed in Lagrangian coordinates, smectics are described in terms of the phase 
variable u, which is Eulerian. To compare with Scherk's surface we shall take d = n here 
and in the following. Because ifi must be single valued, the Burgers vector b must be an 
integer multiple of the layer spacing. The (Lagrangian) surfaces of constant phase h(x, y) 
are defined implicitly through 

h{x,y) - u(x,y,h{x,y)) = kd 7 (6) 

where k is an integer. Thus d x h = d x u/(l — d z u) and similarly for d y h. For low-angle 
grain-boundaries, d z u ~ and the derivatives of h and u are equal. This will not be the 
case for large angle grain-boundaries. With this in mind, we will add together the effect 
of a plane of equally spaced screw dislocations with spacing 

A single screw dislocation is described by u\ — [6/(27r)] tan _1 (y/a;) with b = kd: 
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An array of screw dislocations gives: 
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These sums can be performed explicitly via the Poisson summation formula: 

b sm(2ny/l d ) 



x array '1\ d COS 1 1 ( j! 77./'//,/ ) — ('( )S{ 2 / / ,/ j 

q b sinh(27rx//d) 

y array ^ COSh(2lX X / 1 d) ~ COs(2ny/l d ) ' 



(9) 



Note that the superposition of (^) gives d z u aTTScy = 0, and thus w a rray is independent of 
z. However, w a rray is the displacement due only to the screw dislocations - the smectic 
will relax via a smooth displacement field u smoo th(x, y, z). To find ^smooth we consider 
the boundary condition on the smectic: the strain must vanish at x = ±00. Writing 

1 2 

+ w smoo th, the rotationally invariant strain is u zz = d z u — k (Vm) and the 



u = u 



array 



boundary condition becomes: 



•20,T, - (LKuf = (^-) (1.0) 



or d z u = 1 ± v/l — (b/2ld) 2 - Since u should vanish as Id — > 00 (or b — > 0), we have 



d z u = 1 — y 1 — (b/2ld) 2 . Comparing (^) and (|9p while translating from Eulerian to 
Lagrangian coordinates, we find that Scherk's surface is an array of screw dislocations, 
i.e., a TGB with: 

tan(a/2) = (b/2l d ) [1 - (6/2/ d )]" 1/2 , (11a) 
sin(a/2) = ir/l d , (lib) 
h(x, y) = h(x cos(o;/2), y). (He) 

Equation ( p.la|) implies that sin(a/2) = bj (21 d ), a standard result from the theory of twist- 
grain-boundaries, while ( |116|) shows that b = 2n = 2d, confirming that there are strength 
+2 dislocations. Finally, (|llcj ) displays the only difference between Scherk's surface and the 
screw dislocations: Scherk's surface is dilated along the twist axis by cos(q;/2). We also note 
that the Eulerian coordinate for Scherk's surface is simply u(x, y, z) = u(x cos(a/2), y, z). 

We now consider the free energies of Scherk's surface and an LSD surface with ar- 
bitrary Burgers vector b = kd. This free energy has two contributions. The first is the 
bending energy. In Lagrangian coordinates it can be written as the product of the free en- 
ergy per surface, the density of surfaces per unit length and the total length of the system. 



J dxdy^/l + (d x h) 2 + (d y h) 2 H 2 , (12) 



We use the Helfrich-Cahn energy for each surface, as has been utilized to study lyotropic 
monolayers ||. We find: 

b 2dcos(a/2) 

where H is the mean curvature of the surface. Note that n/d = K\ is the three-dimensional 
splay modulus of the smectic. Scherk's surface is minimal and thus has H = - which is 
precisely why it proves to be a useful ansatz. The LSD surfaces have a non-zero bending 
energy Using the standard non-linear expression for H ||, we find 

kL z F 3 f sinh 2 (x) sm 2 (y) [cosh(x) - cos(y)]" 3/2 

dxdy— — , (13) 



8d cos(a/2) J ' [(} + p) cos h(£) - (1 - T) cos(y)] 5/2 

where V = tan 2 (a/2). Inspection shows that (|T3| ) is a convergent integral, proportional to 
the cross-sectional area of the sample, and that for small angles ~ b 4 LyL z /(dl^). Thus 
for large we see that the interaction energy between two defects scales as Z^~ 5 . This is 
different from the usual exponential interactions that a linear theory with both director 
and displacement fields |7j would predict. Moreover, for fixed angle a = 2 sm~ 1 (b/2ld) 1 
larger values of b are favored over smaller ones. We note that in addition to the mean 
curvature term in ([l^), one could also consider a Gaussian curvature term of the form 
K J dS K where K is the Gaussian curvature. If the elastic constants do not depend on 
the coordinates (x,y), then the integral of the Gaussian curvature is independent of the 
surface geometry. The case of varying elastic constants has been considered elsewhere 0. 

Since smectic and lamellar layers have a preferred spacing, there is an energy cost 
for layer compression. The compression energy F c is most simply written in terms of the 
Eulerian coordinates u: 



F c = | / d *x 



d z u - - (Vu) 2 
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where the nonlinear terms assure complete rotational invariance. The compression energy 
does not vanish for either LSD surfaces or for Scherk's surface. For Scherk's surface F c is 
proportional to the cross-sectional area of the sample in the yz-plane: 



7rcos(a/2) J 



I ^ 2 

D cos(2#) - -D 2 sin 2 (a/2) sin 2 (2#) 
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(15) 



where D 1 = [cosh(27) — cos(20)] ~ \j 2 + 6> 2 ] for 7,S « 1. For an LSD surface this 
compression energy F c is also proportional to the cross-sectional area: 



BL y L z b 4 r /2 

7T 2 (2i d r 



/ dB d<yD 2 cos 2 (29). (16) 
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Note that the leading terms in the expansions of F c and F c in a for a C 1 are identical 
when b = 2d, i.e., for an LSD surface made of +2 dislocations. Because of the relative sign 



between the two terms in ( |i5|) , it is difficult to determine the relative magnitudes of F c 
and F c for the same twist angle a and b = 2d. However, it is easy to see from ([16]) that 
for small angle, strength +1 dislocations lead to a lower compression energy than Scherk's 
surface. 

More importantly, in both cases (LSD and Scherk's surfaces), the integrals diverge for 
small 7 and 9. This means that all of them must be cutoff at some short distance 2ix / A. ~ d 
as has been observed by Kleman ||. Inside the cutoff there is a defect core in which the 
smectic order parameter must vanish. The "core energy" that we find here is actually a 
nonlinear elastic energy, not to be confused with the usual core energy which arises from 
the order parameter vanishing at the defect and from relaxation of director modes. We 
shall use the term "core energy" to refer to the elastic core, keeping in mind that this is 
only part of the total core energy. The integrals (|l5| ) and (|16|) diverge as (Aid) 2 . Thus the 



dominant contribution to the compression energy for l d ^> A 1 ~ d is 

F c -^^(Al d f. (17) 

This has a simple interpretation: since L y /l d is the total number of dislocations, fll7j ) is the 
sum over individual dislocations of elastic core energies eL z with energy per unit length 
e = Bb^A 2 . Unlike the bending energy, this energy strongly favors small Burgers vectors,. 
In addition, we may use (^5[) or fll6| ) to calculate the leading interaction energy between 
defects. One finds that the interaction energy is cutoff dependent through logarithmic 
terms. Therefore, on dimensional grounds the interaction energy per unit length scales as 
Bb 4 l~ 2 \n(Al d ). The cutoff depend ence of the interaction is somewhat unusual: in linear 
elasticity, defect-defect interactions only depend on the defect separation. Because this 
system is non-linear, no such decomposition can be easily made. 

Since both LSD and Scherk's structures necessarily have core regions, the energy of 
the core itself must also be taken into account. A detailed analysis of the core structure 
would be interesting and would allow for explicit energy comparisons between the Scherk 
and LSD structures. In such a calculation, the core size 2% /A itself would be set so as to 
minimize the total energy. This is similar to the analysis of defects in lyotropic lamellar 



systems [10]. Note that there is only one set of parallel core regions, despite the fact 



that there are two, perpendicular sets of helicoid-like structures in Scherk's surface (|T|) 

6 



Once the core regions are present in one set of helicoids, the other set no longer consists 
of perfect helicoids - they are interrupted by the cores of the first set. Thus, one can 
unambiguously identify a single set of dislocations in the LSD and Scherk structures. This 
corroborates our earlier claim that Scherk's surface contains a single set of true dislocations, 
while the second set of (interrupted) helicoids appears only through the geometry of the 
other defects. Presumably, energy considerations will determine which set of helicoids 
become dislocations. At a = tc/2 the energy is degenerate. However, at smaller a the 
preceding discussion suggests that a system with a larger defect separation will have 
a lower energy. This choice of true topological defects is consistent with the traditional 
construction of grain boundaries in the TGB-A phase |4|]. Finally, we note that Scherk's 
surface has a non-zero compression energy. The true equilibrium twist grain boundary 
will adopt a geometry that is an energetic compromise between bending and compression 
deformations. As a result, it will have a non-vanishing mean curvature, and it will have 
a structure that is identical to neither Scherk's surface nor to any LSD surface discussed 
here. It would be interesting to consider a variational ansatz based on Scherk's surface 
with arbitrary, independent dilations of x and y. 

We have shown that Scherk's surface is an anisotropic dilation of a periodic surface 
constructed of a single set of strength 2 screw dislocations. Furthermore since the lamel- 
lar ground state has a preferred layer spacing, layer compression contributes to the free 
energy of the structure. This breaks the dual mapping between the helicoids. It thus 
follows that Scherk's surface is a twist grain boundary composed of a single set of parallel 
screw dislocations and that the geometry of these defects creates a perpendicular set of 
helicoidal structures in the surface. We have also demonstrated that Scherk's surface has 
a higher energy than an LSD structure built of +1 dislocations for small angles. In the 
case of biphasic materials, the +1 dislocations we consider here would be topologically 
forbidden - they would become +1/2 dislocations. In this case the energetics would be a 
competition between Scherk's surface and the +2 LSD. In either case, a detailed analysis 
of the core structure would be required to make an unambiguous prediction of the most 
stable structure at larger angles, whether it is a distorted Scherk's surface, a distorted LSD 
surface or some interpolation between them. 
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Figure Captions 



Fig. 1. Twist-grain-boundary structure for +2 screw dislocations. The horizontal lines 
are the smectic layers at the grain boundary, while the line at ±a/2 are the layers 
from the boundary. The very dark set of horizontal lines follows a single smectic 
layer across grain boundaries. Note that at the grain boundary, the smectic layers 
are dilated by 1/ cos (a/2). 
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Figure 1 



